Coulomb Confinement from the Yang-Mills Vacuum State in 2+1 Dimensions 
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The Coulomb-gauge ghost propagator, and the color-Coulomb potential, are computed in an ensemble of 
configurations derived from our recently proposed Yang-Mills vacuum wavefunctional in 2 + 1 dimensions. 
The results are compared to the corresponding values obtained by standard Monte Carlo simulations in three 
Euclidean dimensions. The agreement is quite striking for the Coulomb-gauge ghost propagator. The color- 
Coulomb potential rises linearly at large distances, but its determination suffers from rather large statistical 
fluctuations, due to configurations with very low values of fio, the lowest eigenvalue of the Coulomb-gauge 
Faddeev-Popov operator. However, if one imposes cuts on the data, effectively leaving out configurations with 
very low fio, the agreement of the potential in both sets of configurations is again satisfactory, although the 
errorbars grow systematically as the cutoff is eliminated. 



I. INTRODUCTION 

Non-perturbative properties of non-abelian gauge theories, 
in particular color confinement, chiral symmetry breaking, 
and the existence of the mass gap, must be encoded in the 
structure of the ground state (vacuum) of these theories. It 
is natural, then, to look for evidence of these properties in 
the vacuum wavefunctional ^o[^] of quantized gauge field 
theory in some physical gauge. The simplest non-trivial set- 
ting is SU(2) gauge theory with no dynamical matter fields, 
and in one lower space dimension where the theory is super- 
renormalizable. 

Recently, we have proposed a very simple approximate 
form for the Yang-Mills vacuum wavefunctional in temporal 
gauge and D = 2 + 1 dimensions [ 1]{3 
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Here T> 2 is the covariant Laplacian in the adjoint represen- 
tation, whose lowest eigenvalue is Arj, m is a constant with 
dimensions of mass proportional to g 2 ~ and B a (x) — 
F± 2 (x) is the color magnetic field strength. On the lattice, 
—T) 2 is given by 
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where the Uk{x) are the link fields in the adjoint representa- 
tion. 



The wavefunctional in Eq. (HJ is reminiscent of 
Samuel's 10], the difference being that in his proposal a 
single free parameter m§ replaces our (— Ao + m 2 ) in the de- 
nominator. The reason for subtracting the lowest eigenvalue 
from the operator (— T> 2 ) is that our numerical simulations 
indicate that the spectrum of this operator may be divergent 
in the continuum limit |Q]] 

The proposed vacuum wavefunctional (Q]) has quite a few 
attractive properties: 

1. In the free-field limit (g — > 0), the covariant Lapla- 
cian turns into an ordinary Laplacian, Ao and m vanish, 
and [A] becomes the well-known vacuum wavefunc- 
tional of electrodynamics: 




d 2 xd 2 y [dxAl(x) -d 2 A a l (x)\ 



[d 1 Al{y)-d 2 A\{y)]\. (3) 



2. The expression (|TJ is a good approximation to the true 
vacuum also in a completely different corner of the con- 
figuration space, namely if we restrict to fields constant 
in space and varying only in time. In the large-volume 
limit the solution of the Yang-Mills Schrodinger equa- 
tion in that case is, up to 1 /V corrections: 



= exp 
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(4) 



and exactly the same expression follows from ([TJ as- 
suming IffAi^l > m, V^o- 
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1 Expressions below are assumed to be properly denned on a lattice, with 
lattice spacing serving as regulator, but for simplicity we will often use 
continuum notations. 



2 It is a little difficult to compare our wavefunctional directly with that of 
Karabali and Nair [3, 4], because their proposed wavefunctional in new 
variables, when converted back to temporal gauge and the usual variables 
of gauge theory, is not gauge invariant. One must therefore suggest some 
gauge invariant extension. The simplest such extension, proposed in 01 
and investigated numerically in j 1], does not seem to give the correct string 
tension. Other extensions are, however, possible. 
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3. If we divide the field strength B a (x) into "fast" and 
"slow" components, the part of the (squared) vacuum 
wavefunctional that depends only on £> s i ow is 



l*r. 
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Such a form is expected on the basis of dimensional- 
reduction arguments JH-0]: it is exactly the probabil- 
ity measure for Yang-Mills theory in two Euclidean di- 
mensions, which (i) is confining for m > 0, and (ii) 
exhibits Casimir scaling for string tensions of all color- 
charge representations. The fundamental string tension 
is easily computed as <jf{P) — 3m/(4/3). The last ex- 
pression can be used to fix the value of m in Eq. ([TJ at 
a given j3 from the known value of erp(/3). 

4. Confinement requires m to be larger than zero. If one 
takes m in the wavefunctional (Q]i as a variational pa- 
rameter and computes (approximately) the expectation 
value of the Yang-Mills Hamiltonian, one finds that a 
non-zero (and finite) value of m is energetically pre- 
ferred. 

5. If we fix the mass m in the wavefunctional to get the 
right string tension op at a given f3, we can test our 
proposal by calculating e.g. the mass gap of the theory. 
We have proposed a recursive procedure for generat- 
ing independent lattice configurations with the proba- 
bility distribution given by the square of the wavefunc- 
tional © (see Ref. CD] and Sec.|III]for details). We call 
two-dimensional lattice configurations obtained in this 
way "recursion lattices". One can compute observables 
with these lattices, and compare the results with corre- 
sponding values obtained from "Monte Carlo lattices", 
i.e. two-dimensional slices of lattices generated in a full 
D = 3 lattice Monte Carlo simulation. It turns out 
that, given the asymptotic string tension as input, the 
mass gap comes out fairly accurately from our wave- 
functional. The discrepancies with the Monte Carlo re- 
sults of Meyer and Teper [8] for the + glueball mass 
are at the level of at most a few (< 6) percent. 

6. The dimensional reduction form © at large distances 
implies an area law fall-off for large Wilson loops, 
and also Casimir scaling of higher-representation Wil- 
son loops. The question is then how Casimir scal- 
ing turns into iV-ality dependence at large distances, 
i.e. how color screening enters in this setting. There 
are indications that terms needed for color screen- 
ing might be contained in our simple wavefunctional 
and would appear as corrections to the dimensional- 
reduction form (©. 

In this article we proceed further in comparing quanti- 
ties derived from our proposed vacuum wavefunctional, with 
those obtained from full Monte Carlo simulations of the 
Yang-Mills theory, by computing the Coulomb-gauge ghost 
propagator and the color-Coulomb potential. The translation 
between temporal gauge and the minimal Coulomb gauge was 



discussed in detail in Ref. [9], but the conclusion (cf. SecUJi is 
simply that the Yang-Mills wavefunctional in Coulomb gauge 
can be obtained by restricting the temporal-gauge wavefunc- 
tional to transverse gauge fields. In Sec.|III]we review the re- 
cursion procedure used to generate lattice configurations with 
probability weighting given by the square of the wavefunc- 
tional ([T]i. Our numerical results are presented in Sec. lIVI first 
for the Coulomb-gauge ghost propagator in Sec. IIV Al and 
then for the color-Coulomb potential in Sec. IIV Bl The latter 
section also contains a discussion of subtleties in the determi- 
nation of the color-Coulomb potential, due to exceptional lat- 
tice configurations. Our conclusions are briefly summarized 
in Sec. [V] 



H. FROM TEMPORAL TO COULOMB GAUGE 

The proposed wavefunctional, Eq. ((T), is formulated as an 
approximate solution of the Yang-Mills Schrodinger equation 
in the temporal gauge. To compute quantities in Coulomb 
gauge we need to find its Coulomb-gauge equivalent. This is 
an easy task, due to the fact that both the temporal and the 
Coulomb gauge are compatible with a Hamiltonian formula- 
tion and a physical transfer matrix. 

In temporal gauge, Ao = 0, in D = 2 + 1 dimensions the 
continuum Hamiltonian has the simple canonical form 



H = 



d 2 x 



(6) 



where the color-electric fields Ef (x) — — i[S/SAf(x)] areop- 
erators canonically conjugate to the spatial components of the 
vector potential. The temporal gauge still possesses invari- 
ance under space-dependent, time-independent gauge trans- 
formations fi(x). The generator of local infinitesimal space- 
dependent gauge transformations, in the absence of external 
color sources, is given by 



G a (x)=-Y d Vf[A]E b i {x), 



(7) 



i=i 



T>i[A] is the covariant derivative. Physical wavefunctionals 
are required to satisfy the Gauss law constraint 

2 r 

G a (xMA] = £ (S ac di + ge abc A\) jj-^W = °> ( 8 > 

i—l 1 

which means that the wavefunctional must be invariant under 
infinitesimal gauge transformations. 

Due to the local gauge invariance of wavefunctionals in the 
temporal gauge, the volume of the gauge group can be ex- 
tracted from inner products of wavefunctionals 

<*i|* 2 > = / [DA] ^l[A]^ 2 [A] (9) 



via the Faddeev-Popov procedure II 1 Ofl . We parametrize con- 
figurations A by A — n A± , where A± is the representative 
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of A in the minimal Coulomb gauge, i.e. A± is transverse, 
V • A± = 0, and belongs to the fundamental modular region 
(FMR), A± G A; fl denotes the gauge transformation that 
brings the configuration A to the minimal Coulomb gauge. 
The Faddeev-Popov formula then gives 

(*i|*a>= / [DA±] debM[A ± ]*l[A ± ]9 2 [A±], (10) 
Ja 

where _M[j4jJ = — V • 2?[Aj_] is the Faddeev-Popov opera- 
tor, symmetric and positive-definite for Aj_ E A. The right 
side of Eq. ( TTOb is the proper expression for the inner product 
in the minimal Coulomb gauge. This means that in the oper- 
ator formalism, the minimal Coulomb gauge is a gauge fixing 
within the temporal gauge of the remnant local gauge invari- 
ance. The wavefunctional in Coulomb gauge is the restriction 
of the wavefunctional in temporal gauge to transverse fields in 
the FMR: 



where 



Coulomb 



[A±] = *[A A 



A± e A. 



(11) 



(The conditions Aq = and V -A = cannot both be imposed 
at all times, but can be imposed at one fixed time, and this is 
all that is required.) 

The vacuum expectation value of an operator Q in Coulomb 
gauge can then be computed from 



(Q) 



^^jy Coulomb 
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Inverting the Faddeev-Popov gauge fixing takes us to 

(Q) = <*o|Q [ a A] |*o>, 



(12) 



(13) 



i.e. we generate configurations following the probability dis- 
tribution ^q, transform them to the Coulomb gauge, and eval- 
uate the observable Q in the transformed configuration. From 
the path-integral representation of the vacuum state, we may 
also go from ( fT3l l to 



(Q) = (Q n 'A(x,t = o) 



(14) 



where the right hand side is the expectation value obtained in 
D = 3 Euclidean dimensions, and f2' is the gauge transfor- 
mation which takes the gauge field on a t = time-slice into 
Coulomb gauge. 



III. GENERATION OF "RECURSION LATTICES" 

To assess the consequences of the proposed wavefunc- 
tional, Eq. (Q~|i, one would like to generate lattice configura- 
tions with probability distribution given by the square of the 
wavefunctional 



P[A] 



exp 



_J_ J d 2 xd 2 yB"(x)1C? y [A}B b (y) 
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We have proposed, in Ref. JH], an iterative procedure to 
achieve this goal. Here we summarize, for completeness and 
reader's convenience, its most important points. 

We define a probability distribution of gauge fields A with 
the kernel JC depending on the background field A 1 : 



P[A;JC[A'}} 



exp 



-i J d 2 xd 2 yB a (x)JCf y [A'}B b (y) 



(17) 



Assuming the variance of the kernel JC is small among ther- 
malized configurations gauge-fixed to an appropriate gauge, 
we can approximate P[A] by 



P[A]^P[A, (/C)] = P 



A, J DA' K{A']P[A'} 



(18) 



DA' P [A, JC[A']] P[A'], 

and solve the equation iteratively 

PW[A] = P[A;1C[0}}, 



P (n+1} [A}= J DA' P[A;K{A']]P^[A'}. (19) 

In practice, we work on a two-dimensional lattice in an ax- 
ial gauge (Ai — 0), which enables one to change variables 
from A\ and A2 to B cheaply, without introducing a field- 
dependent jacobian. Initially, we set also A2 = 0. Then the 
iterative procedure consists of the following steps: 

(i) Given A 2 , set A' 2 = A 2 . 

(ii) The probability P [A; K[A']] is gaussian in B, diago- 
nalize K [A 1 ] and generate new _B-field stochastically. 

(iii) From B, calculate A2 in axial gauge, and compute ev- 
erything of interest. 

(iv) Go back to step (i), repeat as many times as necessary. 

Lattice configurations generated by this procedure are referred 
to as "recursion lattices." It turns out that the procedure con- 
verges rapidly, after O(10) iterations (cycles above), and the 
assumption about a small variance of K, among configurations 
is supported a posteriori by the absence of large fluctuations 
of the spectrum of JC evaluated on individual recursion lat- 
tices. 



IV. COULOMB-GAUGE RESULTS 



3 We have absorbed the coupling g into the definition of A^, which accounts 
for the factor of 1/g 2 in the exponent in Eq. )15t . 



Confinement exists, of course, in any gauge, but in some 
gauges the phenomenon may be easier to understand than in 
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others. Coulomb gauge has received some attention, follow- 
ing the seminal works of Gribov ifTHl and Zwanziger 02 > 
who argued that the low-lying spectrum of the Faddeev- 
Popov operator in Coulomb gauge probes properties of non- 
abelian gauge fields that are crucial for the confinement mech- 
anism. The ghost propagator in Coulomb gauge and the color- 
Coulomb potential are directly related to the inverse of the 
Faddeev-Popov operator, and play a role in various confine- 
ment scenarios. In particular, the color-Coulomb potential 
represents an upper bound on the physical potential between 
a static quark and antiquark, which means that a confining 
color-Coulomb potential is a necessary condition to have a 
confining static quark potential H 1 3fl . It is therefore an im- 
portant check on the validity of our proposed vacuum wave- 
functional, to see how well it can reproduce the values of 
Coulomb-gauge observables, such as the ghost propagator and 
color-Coulomb potential, that can be obtained by standard lat- 
tice Monte Carlo techniques. 



even apart from Gribov copies, the Coulomb-gauge condition 
does not fix the gauge completely. There is a remnant global 
symmetry such that if a set of link matrices Uk(x) satisfies 
the condition, so does U' k {x) = fi Uk(x) Ci>, where fi is 
a space-independent SU(2) group element. The inversion of 
the Faddeev-Popov operator, Eq. ( l22l i. needed to compute the 
ghost propagator d20l l, is therefore performed in the subspace 
orthogonal to the trivial zero eigenmodes. 

At lattice coupling (3 = 6 on a 24 2 lattice, and (3 = 9 
on a 32 2 lattice, we have computed the Coulomb-gauge ghost 
propagator separately in two ensembles of lattice configura- 
tions generated in two different ways: 

(i) recursion lattices are generated by the procedure de- 
scribed in Section [nil me mass parameter m in the 
wavefunctional ((TJ, at each lattice coupling, was cho- 
sen to reproduce the value of the string tension at the 
same coupling as given by Ref. [8|]; 



A. Ghost propagator 

The ghost propagator in Coulomb gauge is given by the 
inverse of the Faddeev-Popov operator: 



g{r)=Um[a\ 




(20) 



x—y\=R 



In D = 2 + 1 dimensions, with SU(2) gauge group and lattice 
link matrices parametrized by 

Un(x) = b^x)! + ia c afa), b^xf +a c fi (x) 2 = 0, (21) 

the lattice Faddeev-Popov operator is given by (assuming the 
lattice version of the Coulomb-gauge condition V • A = is 
satisfied) 

2 
k=l 

- K, v -M*)-s y , x -My)} ( 22 ) 



k=l 



This operator is symmetric and positive-definite anywhere in 
the Coulomb-gauge Gribov region, i.e. for configurations that 
are local minima of the quantity 



(ii) Monte Carlo lattices are generated by Monte Carlo 
simulations of euclidean SU(2) lattice gauge theory in 
D = 3 dimensions with the standard Wilson action. 
From each thermalized configuration only a single (ran- 
dom) space slice at fixed euclidean time was chosen. 
The probability measure for such lattice time-slices, 
when transformed to Coulomb gauge, is given by the 
true Coulomb-gauge vacuum wavefunctional (i.e. the 
vacuum state of the corresponding lattice transfer ma- 
trix). 

There were 1000 lattice configurations in each ensemble, at 
each of the couplings f3 = 6 and 9. 

Each two-dimensional lattice configuration was fixed to the 
Coulomb gauge by minimizing the quantity 7Z, Eq. ( 1231 . via 
the usual (over)relaxation methodQ Then, the inverse of the 
Faddeev-Popov operator was computed using the standard 
linear-algebra tools (Octave). This enabled us to determine 
the ghost propagator directly in the coordinate representation, 
in contrast to most other lattice investigations which generally 
determine propagators in the momentum representation (see 
e.g. Refs. ifilfll). 

The results for our largest coupling (/3 = 9) and lattice size 
(32 2 ) are displayed in Fig. Q] (right). The agreement between 
recursion and Monte Carlo lattices is really quite striking, with 
the differences between the displayed data sets being compa- 
rable to the size of the symbols. The same agreement is ob- 
served also for /3 = 6 on 24 2 lattice, cf. Fig. Q] (left). 



K[U] 



£5>[[/ fc (a;)]. 



(23) 



x fc=l 



More precisely, on a lattice with periodic boundary condi- 
tions the Faddeev-Popov operator possesses three trivial zero 
eigenvalues, the eigenmodes being independent of x. The 
existence of these modes is a consequence of the fact that, 



It should be noted that this procedure returns a gauge copy in the Gribov 
region; no procedure for fixing to the fundamental modular region exists. 
However, we feel that there is nothing sacred (or more physical) about 
Gribov copies in the fundamental modular region, and will be satisfied 
with local, rather than global minima of 1Z. 
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FIG. 1. The Coulomb-gauge ghost propagator: (left) j3 = 6 on 24 2 lattice, (right) j3 — 9 on 32 2 lattice. 



B. Color-Coulomb potential 



The potential between a static quark and antiquark located 
at points x and y, respectively, is proportional to{j 



V(R) 



(M[A]-i(-V 2 )M[A]-% 
1 , _„ 1 



\x—y\=R 



(24) 



y-v[A] 



Provided we know the inverse of the Faddeev-Popov operator 
in a configuration (on the subspace orthogonal to trivial zero 
modes, see the preceding section), the computation of the po- 
tential is quite straightforward. The result is shown in Fig. [2] 
After what we have seen in Section HV Al the figures for 
potentials come as a surprise. In the case of the ghost prop- 
agator, there was almost no difference between recursion and 
Monte Carlo lattice ensembles. Now we observe quite strong 
differences. The origin of these differences can fortunately 
be identified. Both in the Monte Carlo and recursion ensem- 
ble there exist "exceptional" configurations in which the low- 
est nontrivial eigenvalue of the Faddeev-Popov operator jUrj 
is still positive, but extremely small. It means that in the 
space of gauge-equivalent configurations there exist a valley 
along which the minimized quantity 1Z, Eq. d23l ), almost does 
not change. Consequently, these configurations were always 
rather difficult to gauge-fix to Coulomb gaugeQ The existence 
of these exceptional configurations does not have a crucial im- 
pact on the ghost propagator, since its definition contains a 



5 We have omitted below the normalization factor g 2 CV /d,A, where Cf is 
the eigenvalue of the quadratic Casimir operator in the fundamental color 
representation and cLa is the dimension of the adjoint representation of the 
color group, and nevertheless call the above quantity the color-Coulomb 
potential. The normalization factor would be needed if we intended to get 
a numerically accurate value of the Coulomb string tension, which was not 
our goal in the present study. 

6 A similar type of "exceptional" hard-to-gauge-fix configurations was en- 
countered e.g. by the Berlin group 11611 - 



single power of the inverse Faddeev-Popov operator, but their 
influence is strongly amplified in color-Coulomb potentials, 
where the inverse appears twice. 

We will illustrate this point on the data for /3 = 9 on 32 2 
lattice. One can evaluate the values of the potential in a single 
lattice configuration. The lowest eigenvalue of the Faddeev- 
Popov operator directly influences the absolute value of the 
potential at the origin, |V(0)|, so one can classify configura- 
tions by their values of |V(0)|, and evaluate average poten- 
tials from sets of configurations satisfying a number of cuts: 
{|V(0)| <K h i = 1,2,... , K} . 

Figure [3] (left) displays values of |V(0)| in the individual 
configurations from the recursion and Monte Carlo ensem- 
bles. One can see that the majority of configurations lie in 
the range between about 2 to 20, but there are rare instances 
of configurations above 100 (among Monte Carlo lattices) or 
even 400 (among recursion lattices). The frequency of these 
"exceptional" configurations in the ensemble is very small, 
but such configurations give rise to tremendous fluctuations in 
the measured average values of the potential. The distribution 
of configurations according to | V(0) | is further illustrated in 
Fig. fright). 

The K-dependence of the average value of the magnitude of 
the color-Coulomb potential V(R) at R = 0, evaluated in sub- 
sets of Monte Carlo and recursion lattices with |V(0)| < K, 
is displayed in Fig. [4] At j3 = 9 (right panel), the aver- 
age values in both ensembles agree reasonably at least for 
K < O(10). Such a cut on |V(0)| is satisfied by about 85% 
Monte Carlo lattices, and almost 80% recursion lattices (cf. 
Fig. [5J right). The severity of the problem of "exceptional" 
configurations is seen even more dramatically in the data for 
P = 6 on 24 2 lattice (Fig. [4] left panel). Here the close 
agreement of the average | V(0) j between Monte Carlo and 
recursion lattices persists for the cut-off k at least 0(100) 
(satisfied by 99% of lattices); then, a single recursion lattice 
with extremely high value of | V(0)| (i.e. extremely small low- 
est eigenvalue of the Coulomb-gauge Faddeev-Popov opera- 
tor) completely distorts the picture, and causes the disagree- 
ment between color-Coulomb potentials and the huge error- 
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FIG. 2. The color-Coulomb potential computed from all measured configurations: (left) /3 — 6 on 24 2 lattice, (right) /3 — 9 on 32 2 lattice. 
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FIG. 3. Left: |V(0)| in the individual configurations from the recursion and Monte Carlo ensembles. Right: The fraction of configurations 
with | V(0) < k. The top curve is for Monte Carlo lattices, the bottom one for recursion lattices. The area between the curves is colored. 



bars seen in Fig. [2] 

If the color-Coulomb potential is evaluated in Monte Carlo 
and recursion lattices with the same (not too high) cuts ap- 
plied in both ensembles, the fluctuations due to rare configu- 
rations, observed in Fig. [2] are tamed. In Fig.|5]we show the 
results for (3 = 6 with k = 100, and for f3 — 9 with n = 10. 
The potentials are linearly rising over a certain range of dis- 
tances, and agree quite well in recursion and Monte Carlo lat- 
tices. The agreement deteriorates somewhat when the value of 
the cut k is increased, but we believe that approximate agree- 
ment would be restored with sufficient (but obviously huge) 
statistics if the sets of configurations in high |V(0)| bins were 
equally populated. 



V. CONCLUSIONS 

The results presented in this paper strengthen our confi- 
dence that the vacuum wavefunctional for the temporal-gauge 
SU(2) Yang-Mills theory in D = 2 + 1 dimensions, Eq. (Q}, 



is a fairly good approximation to the true ground state of the 
theory. The evidence supplied by Ref . [ 1 ] and summarized 
in Section|T]has been now augmented by measurement of two 
quantities which play a crucial role in understanding confine- 
ment in Coulomb gauge: 

1. The ghost propagator in Coulomb gauge, computed 
from an ensemble of recursion lattices (derived from 
our approximate wavefunctional), and from an ensem- 
ble of Monte Carlo lattices, comes out to be practically 
identical in both ensembles. 

2. If the same cuts on "exceptional" configurations are ap- 
plied in both ensembles, then the color-Coulomb po- 
tential from recursion lattices is also very close to that 
determined from Monte Carlo lattices. Both potentials 
are linearly rising over a range of distances (until the 
effects of lattice periodicity become important). How- 
ever, one would need to considerably increase the statis- 
tics to ensure approximately equal population of excep- 
tional configurations in both ensembles, to convincingly 
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FIG. 4. The average of |V(0)| evaluated in subsets of configurations satisfying the condition |V(0)| < k: (left) /3 — 6 on 24 2 lattice, (right) 
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prove that the deviation between potentials in these two 
ensembles stays tiny even after the cuts are removed. 

While we do not claim that Eq. (Q]i is exact (it is surely only 
an approximation to the true ground state), it has now passed 
a number of very non-trivial tests. More tests are in progress, 
and will be reported on at a later time. 

The extension of our proposed vacuum wavefunctional 
to D = 3 + 1 dimensions is straightforward; it simply 
involves replacing the product B a (x)B b (y) in Eq. (Q~|) by 
Ffj(x)Fjj(y), and replacing two-dimensional integrations by 
three-dimensional integrals. It remains true in 3 + 1 dimen- 
sions that the resulting wavefunctional is exact in the free-field 
limit, and approximately solves the Yang-Mills Schrodinger 
equation in the zero-mode, strong-field limit (see Ref. for 
details). However, going to three space dimensions brings up 



complications associated with the Bianchi constraint. Because 
of that constraint, numerical simulation of the wavefunctional 
along the lines discussed in SectionHIIlis much more challeng- 
ing in D = 3 + 1 than in D = 2 + 1, and a different approach 
to testing our proposal is probably requiredQ 
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7 For some recent work using the dimensionally-reduced form of the 3+1 
dimensional vacuum wavefunctional, cf. Ref. 11711 . 
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